Microscopic calculations of isospin-breaking corrections to superallowed /3-decay 
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The superallowed /3-decay rates that provide stringent constraints on physics beyond the Standard 
Model of particle physics are affected by nuclear structure effects through isospin-breaking correc- 
tions. The self-consistent isospin- and angular-momentum-projected nuclear density functional the- 
ory is used for the first time to compute those corrections for a number of Fermi transitions in nuclei 
from A = 10 to A = 74. The resulting leading element of the CKM matrix, |Kd| = 0.97447(23), 
agrees well with the recent result of Towner and Hardy [Phys. Rev. C 77, 025501 (2008)]. 

PACS numbers: 21.10.Hw, 21.60.Jz, 21.30.Fe, 23.40.Hc 



S 1 



o 

=5 

> 
OS 

m 

On 
O 



X 



Nuclear ft decays provide us with the crucial infor- 
mation about the electroweak force and constraints on 
physics beyond the Standard Model [HQ- Of particular 
importance are superallowed Fermi transitions between 
the J 7r = + members of an isospin multiplet that can be 
used to test the conserved vector current (CVC) hypoth- 
esis and provide the most restrictive test of the unitarity 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Un- 
der assumptions of zero energy transfer and pure isospin, 
the transition matrix elements for superallowed ft decays 
do not depend on nuclear structure. 

For actual nuclei, however, small corrections to the 
Fermi matrix element o fJ = 0+,T = l->J = + ,T=l 
superallowed transitions must be applied (see Refs. 
and Refs. quoted therein): 



\M 



(±)|2 



= 2{l + 5' n )(l + 6 NS -S c ), 



(1) 



where 5c is a nuclear-structure-dependent isospin- 
breaking correction, and 5'^ and 5ns & re radiative correc- 
tions. The corrected product of statistical rate function 
/ and partial half-life t can be written as 
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with Tt being nucleus independent. 

In spite of theoretical uncertainties in evaluation of 
radiative and isospin-breaking corrections, the superal- 
lowed /3-decays provide a stringent test of the CVC hy- 
pothesis. In turn, it is also the most precise source of 
information on the leading element V u a of the CKM ma- 
trix 0] • Indeed, with the CVC hypothesis confirmed, 
V u d can be extracted from the data by averaging over 13 
precisely measured superallowed ft transitions spreading 
over a broad range of nuclei from A=10toyl = 74Q. 

The main focus of this work is isospin-breaking correc- 
tions 5c- This topic has been a subject of numerous the- 
oretical studies using different techniques 

minn. The 



standard in this field has been set by Hardy and Towner 
(HT) 0, H EH who employed the nuclear shell model 
(SM) to account for configuration mixing and the mean- 
field approach to describe the radial mismatch of proton 
and neutron single-particle (s.p.) wave functions. Our 
approach to 5c is based on the self-consistent isospin- 
and angular-momentu m p rojected nuclear density func- 
tional theory (DFT) This framework can si- 
multaneously describe various effects that profoundly im- 
pact matrix elements of the Fermi decay; namely, sym- 
metry breaking, configuration mixing, and long-range 
Coulomb polarization. It should also be noted that our 
method is quantum-mechanically consistent (see discus- 
sion in Ref. [Hj]) and contains no adjustable free param- 
eters. 

The isospin- and angular-momentum projected DFT 
approach is based on self-consistent states \tp) which, in 
general, violate both rotational and isospin symmetries. 
While the rotational invariance is broken spontaneously 
frR Hij . the isospin symmetry is broken both sponta- 
neously (on DFT level) and directly by the Coulomb 
force. Consequently, the theoretical strategy is to restore 
the rotational invariance, remove the spurious isospin 
mixing present in the DFT wave function, and retain 
only the physical isospin mixing caused by the Coulomb 
interaction. This is achieved by the rediagonalization of 
the entire Hamiltonian, consisting the isospin-invariant 
kinetic energy and nuclear interaction (Skyrme) terms, 
and isospin-breaking Coulomb force, in a good-angular- 
momentum and good-isospin basis 
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where Prf T and P[f K stand for the isospin and angular- 
momentum projection operators and Af is the normaliza- 
tion factor. In the current version of the model, nuclear 
isospin-breaking interactions and pairing have been dis- 
regarded. 
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The set of states ([3]) is, in general, overcomplete be- 
cause the K quantum number is not conserved. This dif- 
ficulty is overcome by selecting first the subset of linearly 
independent states (collective space), which is spanned, 
for each I and T, by the natural states \<p; IM\TT Z )^ 



that are eigenstates of the overlap matrix 2fjj . Di- 
agonalization of the Hamiltonian in the collective space 
yields the eigenfunctions: 



|n; <p; IM; T z ) 
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where the index n labels eigenstates in ascending or- 
der according to their energies while I, M, and T z — 
(N — Z)/2 are strictly conserved. By construction, vec- 
tors ((4]) are free from spurious isospin mixing. More- 
over, since projection is applied to self-consistent DFT 
solution, a subtle interplay between the Coulomb po- 
larization (that tends to make the proton and neutron 
wave functions different) and the short-range nuclear at- 
traction (acting in exactly the opposite way) is properly 
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taken into account. As discussed in Refs. [15l. Il6ll2lll22|. 
direct inclusion of monopole polarization effect, which 
is crucial for evaluation of isospin mixing in open-shell 
heavy nuclei, excludes all core-based models thus leaving 
us with essentially one choice: the nuclear DFT. Recent 
experimental data on isospin impurities deduced in 80 Zr 
from the giant dipole resonance 7-decay studies [23| are 
consistent with the magnitude of isospin mixing calcu- 
lated with isospin-projected DFT [22j, and this is very 
encouraging. 

As demonstrated in Ref. [21|, in odd-odd N = Z nu- 
clei, the isospin projection alone is not sufficient and a 
simultaneous angular-momentum projection is a must. 
Unfortunately, this leads to the appearance of singular- 
ities in the energy kernels (22[, thus preventing us from 
using modern Skyrme energy density functionals (EDFs) 
as none of them is usable, whereas those depending on 
integer powers of the density, which are regularizable 
2J], are not yet developed. Hence, at present, the only 
practical option is to use the Hamiltonian-driven EDFs 
which, for Skyrme-type functionals, leaves only one op- 
tion: the density-independent SV parametrization [25] 
supplemented by tensor terms. 

The unusual form of SV impacts negatively its overall 
spectroscopic quality by impairing such key properties as 
the symmetry energy [22j, level density, and level order- 
ing. These deficiencies affect the calculated isospin mix- 
ing. For instance, for the case of 80 Zr discussed above, 
SV yields the isospin mixing 2.8%, i.e., smaller than the 
mean isospin mixing 4.4% averaged over nine commonly 
used Skyrme EDFs, see Fig. 1 of Rcf. 21]. Of course, for 
the description of 5c, of importance is not the absolute 
magnitude of isospin mixing but its difference between 
parent and daughter states 13J. The lack of reason- 



able EDF is, admittedly, the weakest point of our current 



calculations; nevertheless, no significant improvement of 
this aspect can be expected in the near future. 

The + — > + Fermi /3-decay proceeds between the 
I J = 0,T » 1,T Z = ±1) ground state (g.s.) of the even- 
even nucleus and its isospin- analogue partner \I = 0, T ~ 
1,T Z = 0) in the N = Z odd-odd nucleus. While the 
DFT state representing the even-even nucleus is unam- 
biguously defined, the DFT state used to compute the 
N = Z wave function is the so-called anti-aligned con- 
figuration \<p) = \v ® 7r) (or \v ® 7f}), selected by plac- 
ing the odd neutron and the odd proton in the lowest 
available time-reversed (or signature-reversed) s.p. or- 
bits. The anti-aligned configurations manifestly break 
the isospin symmetry but they provide a way to reach 
the \T w 1,1 = 0) states in odd-odd N = Z nuclei 
(ljj ]. This situation creates additional technical prob- 
lems. The anti-aligned configurations appear to be very 
difficult to converge in the symmetry-unrestricted DFT 
calculations. This can be traced back to time-odd com- 
ponents of the EDF. In fact, only in a few cases were we 
able to obtain symmetry-unrestricted self-consistent solu- 
tions. This forced us to impose the signature-symmetry 
on other DFT wave functions, which implied a specific 
s.p. angular- momentum alignment pattern (26j . 

The calculations presented here were done using the 
DFT solver HFODD (v2.48q) [3], which includes both 
the angular-momentum and isospin projection. The cal- 
culated values of 5c depend on the basis size. In order 
to obtain converged result for 5c with respect to basis 
truncation, we use 10 oscillator shells for A < 40 nuclei, 
12 oscillator shells for 40 < A < 62 nuclei, and 14 oscil- 
lator shells for A > 62 nuclei. The resulting systematic 
errors due the basis cut-off do not exceed ~10%. 

The equilibrium quadrupole deformations (/?2 , 7) of the 
anti-aligned configurations in odd-odd nuclei are, in most 
cases, very close to those obtained for even-even isobaric 
analogs. Typical differences do not exceed A/?2 ~ 0.005 
and A7 « 1° except for nearly spherical systems A = 14 
and A = A2, where the concept of static deformation is 
ill-defined, and for A = 10 and A = 18 pairs where odd- 
odd and even-even partners have fairly different shapes. 
As we shall sec below, such deformation difference results 
in large values of 5c- 

All studied odd-odd nuclei, except for A = 14, 38, 
and 42, are deformed; thus, to carry out projections, we 
could use for them the unique lowest anti-aligned DFT 
states. Also for A = 14 and 38, unique configurations 
based on the lpi/2 and 2si/ 2 subshells were used. A dif- 
ferent approach was used to compute 5c in near-spherical 
A = 42 nuclei. In 42 Sc, four possible anti-aligned DFT 
configurations built on the s.p. orbits originating from the 
spherical I/7/2 subshells can be formed, and the corre- 
sponding DFT states differ slightly due to configuration- 
dependent polarizations 21]. Consequently, to evaluate 
5c for A = 42 we took an arithmetic mean over the values 
calculated for all anti-aligned configurations. 
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TABLE I: Experimental ft- values (in sec); Sc values adopted 
in this work (in %); calculated J-t- values (in sec); empirical 
corrections ([5]) (in %), and individual contributions to \ 2 used 
in the CL test. 

Parent ft Sc Tt ^ EXP) X i 



T, = -1 : 



io c 


3041.7(43) 0.559(56) 3064.8(48) 0.39(14) 


1.3 


U Q 


3042.3(11) 0.303(30) 3072.3(21) 0.38(06) 


1.5 


22 Mg 


3052.0(70) 0.243(24) 3082.2(71) 0.64(23) 


3.0 


34 Ar 


3052.7(82) 0.865(87) 3063.5(87) 0.65(27) 


0.6 


T z = : 

26 AJ 


3036.9(09) 0.494(49) 3066.7(20) 0.39(04) 


6.8 


34 C1 


3049.4(11) 0.679(68) 3069.8(26) 0.67(05) 


0.0 


42 Sc 


3047.6(12) 0.767(77) 3069.2(31) 0.74(06) 


0.1 


46y 


3049.5(08) 0.759(76) 3069.0(30) 0.73(06) 


0.3 


50 Mn 


3048.4(07) 0.740(74) 3068.3(31) 0.69(07) 


0.7 


54 Co 


3050.8(10) 0.671(67) 3073.0(32) 0.77(08) 


1.5 


62 Ga 


3074.1(11) 0.925(93) 3088.7(41) 1.52(09) 


41.0 


74 Rb 


3084.9(77) 2.06(21) 3064(11) 1.88(27) 


0.4 



FIG. 1: (Color online) Calculated isospin-mixing corrections 
Sc for T z = -1 -> T z = (top) and T z = -> T z = 1 
(bottom). Our adopted values (stars) arc compared to values 
of Ref. 15] (dots; including errors) and Ref. [Tl| (triangles) . 



The unusually large correction 6c ~ 10% has been 
calculated for A = 38 nuclei. Most likely, this is a con- 
sequence of incorrect shell structure predicted with SV. 
Specifically, as a result of incorrect balance between the 
spin-orbit and tensor terms in SV, the 2s\/2 subshell 
is shifted up to the Fermi surface. This state is more 
sensitive to time-odd polarizations than other s.p. states 
around 40 Ca core, see Table I in Ref. [27J. Consequently, 
the 38 K— » 38 Ar transition has been excluded from our cal- 
culation of Vud- 

The adopted values of Sc are shown in Fig. [1] for the 
T z = — 1 — > T z = transitions in light systems (upper 
panel) and for T z = — > T z = 1 transitions pertain- 
ing to heavier systems (lower panel). It is instructive to 
compare our results to those of Refs. HI and [ll|. In 
light nuclei, the calculated Sc are sensitive to the local 
shell structure. Indeed, although our values of Sc show 



roughly the same trend as those of Ref. [B| , the individual 
values differ. The reason can be traced back to the poor 
spectroscopic quality of SV, which manifests itself much 
stronger in light than in heavier nuclei due to the low 
s.p. level density. Let us also remind that the equilibrium 
minima in parent and daughter nuclei with A = 10 and 
A = 18 differ, and this results in increased Sc values. As 
verified by DFT calculations using other EDFs, and also 
findings of Ref. |Tl| , with higher level density in heavier 
nuclei the detailed shell structure seems to play a lesser 
role. This indicates that gross features of configuration 



mixing in heavier nuclei associated with long-range time- 
even (shape) correlations are less dependent on a EDF 
parametrization and may be relatively well captured by 
SV. The calculated values of Sc for heavier nuclei are in- 
deed quite consistent with the HT results [f|, with the 
exception of A — 62 and 66. 

The predicted isospin-breaking corrections are listed 
in Table HI All other ingredients needed to compute Tt- 
values from Eq. ([2]), including empirical ft- values and 
radiative corrections 6^ and 6ns, were taken from the 



most recent compilation [28[. In the error budget of 



Tt in Table |H apart from errors of ft and radiative 
corrections, we include 10% systematic uncertainty in 
the calculated Sc due to basis truncation. The aver- 
age value J-t — 3070.4(9) s was obtained using Gaussian- 
distribution- weighted formula to conform with standards 
set by HT. This leads to |Kd| = 0.97447(23) which co- 
incides with both the HT result \V 



(HT) 
ud 



0.97418(26) 

Q and a central value obtained from the neutron de- 
cay \V$\ = 0.9746(19) 0. Combining the calcu- 
lated | V ud | with the values of \V US \ = 0.2252(9) and 
1Kb | = 0.00389(44) provided in Ref. 0, we obtain 
|V ud | 2 + |V US | 2 + |V ub | 2 = 1.00031(61), which implies that 
unitarity of the CKM matrix is satisfied with precision 
of 0.1%! 

While our value of | V u d | is consistent with both HT and 
neutron-decay results, a question arises about its confi- 
dence level, especially in light of poor spectroscopic prop- 
erties of S V. To this end, we carry out the confidence-level 
(CL) test proposed recently in Ref. [28| using variant in- 
cluding uncertainties on experiment, 6^, and 6ns- The 
test is based on the assumption that the CVC hypothesis 
is valid to at least ±0.03% precision, implying that a set 
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FIG. 2: (Color online) Calculated (black dots) and empirical 
(white dots, with error bars) values of So as function of pro- 
ton number in the daughter nuclei. Vertical arrow marks the 
values for A = 62. See text for details. 



of structure-dependent corrections should produce a sta- 
tistically consistent set of Ft values. Since only one set of 
calculated 6ns corrections exists [3j, "empirical" isospin- 
symmetry-breaking corrections can thus be denned by 



.(EXP) 

°c 



1 



Ft 



ft(i+s' R y 



(5) 



and they are tabulated in Table fl] The CL can be as- 
sessed by minimizing the root mean square deviation 
between predicted and empirical values of 6c with re- 
spect to Ft in Eq. ([5]). The final result corresponding to 
Ft = 3070.0 s is shown in Fig. [51 Individual contributions 
to x 2 are also displayed in Table Q] Our value of reduced 
X 2 (per degree of freedom; in our case rid = LI) is 5.2. 
This is considerably higher than the values reported in 
Ref. [Il for the Damgaard model [lEl (1.7), SM with 
Woods-Saxon radial wave functions [5| (0.4), SM with 
Hartree-Fock (HF) radial wave functions [4j,l9| (2.2), and 
relativistic Hartree plus RPA model of [LL| (2.1). The 
low CL of our model results primarily from the single 
point at A = 62. 

In summary, the state-of-the-art isospin- and angular- 
momentum-projected DFT calculations have been per- 
formed to compute the isospin-breaking corrections to 
+ — > + Fermi superallowed /3-decays. Our results for 
Ft = 3070.4(9) s and |Kd| = 0.97447(23) were found to 
be consistent with the recent HT value [5] . While the CL 
of our Sc values is low, primarily due to a poor spectro- 
scopic quality of the EDF used, our framework contains 
no adjustable parameters and is capable of describing mi- 
croscopically all elements of physics impacting 6c- The 
results presented in this paper should thus be considered 
as a microscopic benchmark relative to which the further 
improvements (e.g., regularizable EDF and/or inclusion 



of pairing) will be assessed. 
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